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The integral equation of Hallen’s type is derived in the time domain. The solution of the 
integral equation is carried out numerically, and the electromagnetic behavior of linear 
antennas and scatterers under various excitations is presented. The integral equation may 
be used to obtain time-domain responses of coupled parallel linear antennas, and scatterers 
with loads. Combining with the method of characteristics in solving transmission line prob- 
lems, the integral-equation formulation is applied to the cases in which antennas and scat- 


terers are connected to transmission lines. 


INTRODUCTION 


Recent development in computational techniques 
has prompted an increasing interest in the direct 
time-domain investigation of the transient responses 
of antennas and scatterers, Indeed the time-domain 
formulation of electromagnetic problems is more 
fundamental and often offers more physical insight 
than the familiar steady state approach, and some- 
times it also provides computational advantages. 
Transient behavior of linear antennas has been in- 
vestigated in terms of superposition of steady state 
solutions [Schmitt et al., 1966; Tesche, 1971] and 
by space-time-domain integro-differential equations 
[Bennett and Auckenthaler, 1971; Sayre and Har- 
rington, 1972; Miller et al., 1973]. 

Quite recently, a new technique called the singu- 
larity expansion method [Tesche 1973] has been in- 
troduced. This method describes the properties of the 
wire by the location of the singularities of the re- 
sponse function in the complex frequency plane. 
Once these locations and the corresponding residues 
are known, the time responses for any change in the 
angle of incidence or the time behavior of the inci- 
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dent wave can be easily computed. In this paper, we 
shall present an alternative approach which results 
in a Hallen’s type space-time domain integral equa- 
tion; i.e., there are no space-time derivatives of the 
unknown in the integral operator. As will be out- 
lined in the following paragraph, this method has 
definite advantages over other methods. 

The superposition method is time-consuming com- 
putationally: for a simple antenna, it takes over 400 
sec on a CDC 6600 computer to calculate the time 
response of the current on the antenna; whereas the 
present method takes less than 10 sec on a CDC 
6400 computer, which is a considerably slower ma- 
chine. The superposition method also suffers inac- 
curacy due to the convergence problem of super- 
imposing many solutions. The singularity expansion 
method is not very satisfactory in the evaluation of 
the very early time behavior, and typically it takes 
about 250 sec on a CDC 6600 computer to search 
for the singularities and construct the solution. In the 
absence of information on computational time of 
the space-time domain integro-differential equation 
methods, it suffices to point out that the integro- 
differential equations contain time and/or space de- 
rivatives, and the numerical representations of these 
derivatives are not satisfactory; these methods, as 
presented, involve complicated numerical procedures 
and hence are computationally less efficient. Some- 
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times, the solutions contain wiggles [Sayre and 
Harrington, 1971], indicating an unsatisfactory nu- 
merical scheme. Due to the absence of space-time 
derivatives, the numerical procedure in solving the 
integral equation is considerably faster than those 
offered by previous investigators, and we can then 
afford to compute complicated configurations of 
antenna arrays. 

In a recent paper, Poggio [1971] has presented 
the space-time domain integral equation by the 
Fourier transform of the familiar frequency domain 
integral equation. We shall present the derivation of 
the time-domain equation for linear antennas di- 
rectly from the space-time potential functions, which 
appears to be very fundamental in the integral- 
equation formulation of Maxwell’s equations. The 
mathematical steps of deriving the space-time in- 
tegral equation also adds new physical insight to the 
understanding of linear antennas. 

The numerical procedures of solving the space- 
time integral equation is presented and some inter- 
esting results are given, All indications show that the 
space-time integral-equation approach offers accurate 
results and some computational advantages. The 
method is then extended to solve for problems in- 
volving coupled parallel antennas. Finally, we com- 
bine this method with the method of characteristics 
in solving transmission-line problems. by Branin 
[1967] and Liu [1968] to deal with the case of a 
thin wire which is driven by means of a transmission 
line and the problem of parallel antennas which are 
connected by transmission lines such as in the case 
of log-periodic dipole antennas. 


MATHEMATICAL FORMULATION 


The derivation of the space-time domain integral 
equation proceeds in parallel to that of Hallen’s 
integral equation in the space-frequency domain. 
Let us consider a z-oriented wire and start with the 
wave equation [Stratton, 1941] relating the z com- 
ponent of the magnetic vector potential A,(z, ¢) and 
the z-directed incident electric field Ej,.(z, t), 


[0° A,(z, )]/dz’ — c?[8? A(z, OY/ar 
= —l0E;,.(z, t)]/dt (1) 


where c and « are the velocity of light and the dielec- 
tric constant of the medium, respectively. Equation 
1 is valid on the surface of a perfect conductor. 

A particular solution of equation 1 can be ob- 
tained by using the new variables, 


n=ct—z 
E=ct+z 
In terms of é and y, equation 1 becomes, 
[07/dE 9n)14,€, 1) = (1/4)e 
[OE ine)/9tE — 0)/2, E + 0)/2] (2) 


The particular solution of equation 2 is, 


4 § 
aren = 1/4 ff 
[@Eine)/Ot1IG’ — 0/2, @ + 1/21 dt! da’ 3) 


Changing the variables back to the (z, ¢) coordinates 
we get 


Art, 1) = 1/4) [f  lloEinte’, Var") 


- {19@’, 9’)/[8@’, #)]} dt’ de’ (4) 
where the Jacobian is given by 


= 2c 


1 -1 
e o¢ 


[a¢’, n’)1/19@’, t’)] = 


The limit of integration in equation 4 should now 
be changed to (z, ¢); this procedure is illustrated 
in Figure 1. The upper limits in equation 3 give 


Sa 3 
f= 7 
which in the (z, t) coordinates become 
ce! +2’ =ct+z 
ce! — 2’ =ct—2z 
If we integrate / first in equation 5 the two limits 
are given by the straight lines 


Y=rt+z-—2je 
f=t-(Zz-—2/e 
The integration in equation 4 now becomes 


aren=@/yf ff —— SEiol TED ay ae! 
+/af f ek {[ Binc(2’, #”)/a0"} dt! de’ 
= (@/2 J Eindle’, t— @ — 2c] ee! 
+ (c/2) [ Ejnolz’, t + (& — 2’)/e] da’ 


= (ec/2) [ Einc(z’, t — |z — 2’|/c) dz’ (5) 


t'=t -(z-z)/c t'=t+(z-z')/c 


z-ct 


Fig. 1. The limits of integration for A,’ in (z’, ft’) space. 
Since the differential equation (1) is only limited 
to the surface of the wire, we have 


L 
A? @, 1) = (2Z)" i. Ejno(2’, t — |z — z’|/c) de’ (6) 


where L is the total length of the wire, and Z) = 
(ec) the intrinsic impedance of the medium. 
Taking the homogeneous solutions to the wave 
equation (1), and expressing the vector potential 
in terms of the z-directed surface-current density 
J,(z, t) along the wire [Stratton, 1941], we have 


[ff tae, t= be = Poe Ie — Dp) as" 
iE 
= (22 [ Ej(z!, t ~ |z — 2'|/c) de! 


+ filet — z) + fact + 2) (7) 


where Sis the surface supporting J,(r, f), and f,(ct — z) 
and f,(ct + z) are the homogeneous solutions deter- 
mined by the boundary conditions that J,(0, 4) = 
J(L, ) = 0. 

When the observation and source points are far 
apart, it is possible to use the thin-wire approxima- 
tion for the kernel in equation 7, with I(z, t) = 
2a J,(z, t), where a is the radius of the thin wire, 
and we obtain 


E 
[ UG’, t — fe — 2"\/c)V/{4el@ — 2! + PY} ce! 


L 
= (2Zo)* i; Eino(2’, t— |z = z’\/c) dz! 
Oo 


+ filet — z) + falet + z) (8) 
Either equation 7 or 8 can be used for the solu- 
tion on the thin-wire problems by means of nu- 
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merical techniques on a digital computer. The 
physical interpretation of either equation 7 or 8 
helps directly towards the procedures of numerical 
solutions and will be discussed more fully. 


THE NUMERICAL METHOD AND RESULTS 
Equation 7 or its thin wire version, equation 8, 
indicates that for a specific pair of values (z, to), 
the current I(z’, to — |z — z’|/c) lies on the straight 
lines 


t-—@-—27)/c 


t+@—2')/c 
in the (z, t) plane as shown in Figure 2. For a 
discrete set of parameters tp, these straight lines form 
two families of trajectories, 2 and 8, known as the 
characteristic curves [Courant and Hilbert, 1962]. 
On the « family of curves, 


h=t+(z—-zV/e 


to 


to 


or 


de’ /dt =c¢ (9) 
we find current waves propagating along the wire 
in the +z direction. Similarly on the f family of 
curves we find current waves in the —z direction, as 


dz'/dt = —c (10) 
The solution of the homogeneous wave equation 


fi(ct — z) is invariant along an a curve, which is 
characterized by a constant value of ct — z. Similarly 
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Fig. 2. Diagram of z — t showing the two families of 
a and f characteristic curves. Here, Az = L/16. 
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Fig. 3. Driving-point current response for a center-fed 


antenna under a unit-time step voltage excitation of @ = 2 
log L/a = 10, where r = L/C and Y¢ is the retardation time. 


fe(ct — z) is invariant along a 8 curve. These two 
functions should be included in the equations to 
account for the multiple reflection of the current 
at the two ends of the wire. Mathematically, the 
homogeneous solutions enforce the boundary con- 
ditions 1(0, t) = I(L, t) = 0. In effect they sum- 
marize the history of the current by the boundary 
conditions so that the integration may terminate at 
z= Oand2? =L. 

The current at time t may be found from equation 
8 if the current and /,, fe along the e@ and 6 char- 
acteristic curves are known prior to time t. To 
facilitate the determination of the functions f, and fe, 
we extend the characteristic curves beyond the t = 0 
line as shown in Figure 2. Realizing that 


fAct-z=f(t+tz=0 #50 

and that f, and f. are invariant along these char- 
acteristic curves, we find that these functions vanish 
as long as the (z, t) pair is on a characteristic curve 
which extends into the region t < 0. To find the 
current at f = At, we shall first assume the wire to 
be unexcited, i.e., J(z, 0) = 0, and by the previous 
reasoning, we find f,,2(cAt = z) = 0. Thus, the 
current at t = At is readily found from equation 8. 
At t = 2At, since 1(0, t) = I(L, t) = O we find 
fi(2cAt) and fe(2cAt + L) on the outgoing char- 
acteristic curves at z = 0 and z = L, respectively. 
Using similar steps, we may find I(z, t) at t = 2At, 
3At, 4At, --- , and fi(ct) and fe(ct + L) att = 
4At, 6At, ---. The computation is a step-by-step 
time-marching process, based on the condition that 
the currents prior to the point under evaluation are 
known. Hence, if we can perform the calculation 
for one time line, we are able to compute for the 
current at all subsequent time. 


The method of the numerical approximation of 
the integrals in equation 7 or 8 is well-documented 
[Harrington, 1968] and will not be detailed here. 
To achieve better accuracy, the current is inter- 
polated. As observed from the equations, the cur- 
rent, propagating in the directions of the char- 
acteristic curves, should have smooth values along 
them; hence the interpolation is carried out linearly 
along the characteristic curves. It is also worth 
noting that the high frequency components of the 
solution are considerably improved if the singular 
integral, i.e., that at the singular cell 2’ = z, is 
evaluated by the exact current integral with the aid 
of the method of auxiliary integral [Lin and Mei, 
1968]. The difference between the so-called exact 
value and the thin-wire approximation can be more 
than 10% [Liu, 1972]. For nonsingular cells, the 
thin wire approximation is extremely good for 
Az > Sa. 

Some interesting transient responses are pre- 
sented here. In Figure 3, the driving-point current 
responses of a center-fed antenna excited by a unit 
time step voltage is shown for the case Q = 2 log 
L/a = 10. The reflection effects from the two ends, 
the rapid decay of high-frequency components, and 
the eventual no-current flow at late time are clearly 
shown. In Figure 4, the normalized far-field re- 
sponses of the above antenna are presented for 
6 = 15° to 90° at 15° intervals. The far-field formula 


Fig. 4. Normalized far-field responses for a center-fed 
antenna under a unit-time step voltage excitation of Q = 2 
log L/a = 10 where +r — L/c and f is the retardation time. 


is readily derived [Liu, 1972] and is given by 
E,(r, t) = [usin 6)/(4ar)\(0/ar) 
L, 
“ [ Tle’, t- + (2’/c) cos 6] dz’ = (I) 


where ?, is the retarded time to position r. 

When the wire acts as a scatterer with a lumped 
load at z = %, the electric-field integral in equation 
8 has to be modified by introducing an additional 
term, —Vicaa(t — |Z — Z|/c), to account for the 
voltage drop across the load. If the load is an in- 
ductance or a capacitance, Vi..a(t) contains the time 
derivative or the time integral of the load current, 
respectively. The problem is then solved by nu- 
merical approximation of the derivative or the in- 
tegral. In Figure 5, the load currents for a center- 
loaded scatterer under the excitation of a 30° 
incident plane wave are presented. The resistive 
load reduces the current amplitude of the unloaded 
wire. The inductive load increases the time between 
two subsequent zero crossings of the current, in- 
dicating an effective lengthening of the wire; the 
waveform is also smoothed. The capacitive load 
shows the opposite effect and effectively shortens 
the wire; the waveform is sharpened. 

It is worth mentioning that by Fourier trans- 
forming the transient driving-point current and the 
applied time-varying voltage, the input impedance 
of the antenna can be obtained as a function of fre- 
quency. The results so obtained agree very well 
with direct frequency-domain solutions. The time- 
domain approach, however, has a great saving in 
computation time over the frequency-domain ap- 
proach [Liu, 1972]. 


COUPLED PARALLEL ANTENNAS 


The above integral equation formulation can 
easily be extended to the case of coupled parallel 
antennas, For the case of N antennas, we have 


[thee — eZ Vou tate — 2 $a) a 
i 
a > f If’, t —Ie@ — 2’ +.5,7)7/c Ze 
at 


4n[(z — 7) + 6,,°7" 


= (2Z)* i Eine(z’, t — |z — 2’|/c) dz’ 
Lj 


+ fist — 2) + filet tz) f= l,es, N (12) 


where the subscript j.denotes the wire on which z 
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Fig. 5. Induced-load current responses on a center-loaded 
scatterer under excitation of a 30° incident unit-step electric 
field: (2) no load and 50 @ resistive load and (b) in- 
ductive load of j 63 9 and capacitive load of —j 53 Q at 
frequency equals c/L. Here, 9 = 10, and r = L/c. 


is located, by is the perpendicular separation be- 
tween wires i and j, and a, is the radius of wire j. 

Equation 12 is similar to equation 8 except for 
the integrals under the summation sign which ac- 
count for the coupling effect of the neighboring 
antennas. The numerical solution procedure follows 
closely that of a single wire. At time ¢ and at each 
point z on wire j, there is only one unknown, i.e., 
I;(z, t), because the extra integrals contain quan- 
tities that are time retarded due to the separation 
between the antennas i and j, and are known for 
this step-by-step time-marching process. 

To facilitate comparison with the frequency- 
domain results of Chang and King [1968] a pair 
of nonstaggered, identical parallel antennas is 
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SELF {(mmhos) 


Fig. 6. Self-admittance of a pair of nonstaggered, identical, 

coupled parallel antennas separated by half the antenna 

length. Here 9 = 8.5. Legend: conductance; —~— 
susceptance; @ Chang and King [1968]. 


chosen. By applying a voltage to antenna 1 while 
short-circuiting antenna 2 at the driving point, self- 
and mutual-admittances can be obtained by Fourier 
transforming the transient currents and voltage. In 
Figures 6 and 7 the self- and mutual-admittances of 
two antennas with Q = 8.5, separated at a distance 
equal to half the length of the. antennas are pre- 
sented as a function of frequency. More extensive 
comparison is given by Liw[1972]. 


WIRES CONNECTED WITH TRANSMISSION LINES 


In practice, most antennas are driven by means 
of transmission lines, We present here a numerical 


MUTUAL (mmhos) 


Fig. 7. Mutual admittance of a pair of nonstaggered, 

identical, coupled parallel antennas separated by half the 

antenna length. Here Q = 8.5. Legend: conductance; 
-—-— susceptance; ® Chang and King [1968]. 
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Fig. 8. Diagram of x — ¢ plane showing the two families 


of a; and §; characteristic curve. 


procedure in solving such a problem. Essentially, 
the method combines two techniques: the numerical 
method we just outlined for the thin wire and the 
method of characteristics in solving transmission- 
line problems as detailed by Branin [1967] and 
Liu [1968]. 

In summary, the method of characteristics enables 
the transmission line at any point to be represented 
by two Thevenin-equivalent circuits with time-vary- 
ing sources. Each equivalent circuit is defined along 
one set of characteristic curves (a; and f; as shown 
in Figure 8). For a lossless transmission line, the 
equivalent voltage sources, ; 


Sax, t) = o,(x, t) + Roi, 1) 

Sa(x, t) = v.(x, 1) — Roi x, 1) 
are invariant along the respective’ characteristic 
curves; this property is utilized for the numerical 
method. For a lossy transmission line, the method 
is still applicable but is more complicated. 

It is important to point out that this method is 
again a step-by-step time-marching technique and 
is very compatible with the numerical scheme for 
the thin wires. , 


Fig. 9. Diagram of Thevenin-equiv- 

alent circuits for x = 0 at section 

A-A in Figure 8, Ss(x, t) = v:(x, t) 
— Riii(x, t). 


Fig. 10. Diagram of Thevenin-equiv- 

alent circuits for x = d at section 

B-B in Figure 8. Se(x, t) = v(x, t) 
— Riis (x, t). 


When a transmission line is terminated at x = a 
by a thin wire at the position z = 2%, as shown in 
Figure 8, the problem, by using the two methods 
just outlined, reduces to a time-marching solution 
of the two separate networks in Figures 9 and 10. 
In this case, we choose equal time steps for the 
two numerical methods. In Figure 10, the thin wire 
current at z = 2 being the same as i;(d, 1), can 
be obtained by solving the following equation which 
is a modification of equation 8 to take into account 
the Thevenin equivalent circuit, 


L . 
i] LG", t — jz — 2! |/e)/{4el@ — 2/9 + aT} de’ 
L 
. eza| f Eino(e’, ¢ — [2 — 2'\/c) de! 
0 


+ S.(d, t— lz i Za|/c) —_ Rol@o, i= lz ~al/0) | 


+ fi(ct — z) + fa(ct + z) (13) 
The solution procedure is similar to that outlined 
previously. If it is not a scattering case, the integral 
on the right-hand side is not required, In Figure 11, 
we show the result of short-circuit and open-circuit 
termination at x = 0. The scatterer, having Q = 10, 
is excited by a unit time step electric field incident 
at 30° to the z axis. The transmission line has a 
characteristic impedance of 50 © and a transit time 
7; being half of that of the scatterer; it is connected 
to the midpoint of the wire. We observe the pre- 
dominantly capacitive behavior for the short-circuit 
case and inductive behavior for the open-circuit case. 

When parallel wires are connected by transmis- 
sion lines, similar technique is applied. We consider 
here the case of a pair of identical, nonstaggered 
parallel antennas with © = 10, separated by the 
length of the antenna, and connected by a 50-0 
transmission line with transit time 7;. = 7 and driven 
via another 50-0 transmission line with transit time 
t1 = 0.75 +. In Figure 12, we present the re- 
sponses when the system is driven by a unit time- 
step voltage source with a 50-0 source resistance. 
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We observe in Figure 12 that at large t, when 
the antenna currents die down, the transmission-line 
driving-point voltage settles to that of the voltage 
source, because the load to this transmission line 
now behaves open-circuited. 

The frequency domain results obtained by Fourier 
transforming the time responses obtained by the 
above technique compare favorably with those of 
direct frequency domain methods for both single- 
wire cases and multiple-wire cases [Liu, 1972]. 


CONCLUSIONS 


The formulation of the time-domain Hallen’s 
integral equation is of fundamental importance in 
the understanding of the transient behavior of linear 
antennas and scatterers, The equation describes how 
the current at a space-time coordinate (z, t) de- 
pends on previous current values along the wire 
and the incident electric field. It is pointed out that 
the current I(z, t) is directly affected by the cur- 
rent that lies on the two characteristic curves passing 
through (z, t). as well as the previous multiple 
reflections connected with these characteristic curves. 
The condition of the reflections is included in the 
equation by the solutions to the homogeneous wave 
equation, namely, f,(ct — z) and fo(ct + z). 


ANTENNA LOAD-POINT CURRENT {mamp) 


Fig. 11. 


Load point current of a scatterer with Q = 10, 
center-connected to a 50: transmission line (a) with a 
short-circuit load and (b) with an open-circuit load. Transit 
time = 7/2, where r = L/c. 
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ANTENNA DRIVING-POINT CURRENT (mamp) 


VOLTAGE (v) 


8 i/t 
(c)} 


Fig. 12. Antenna driving-point current responses of (a) 

antenna 1 and (b) antenna 2. (c) The transmission-line 

driving-point voltage response. For the whole figure Q = 10, 
7 = L/c, t= 0.75 7, and r1s = 1. 


The numerical method is developed with the idea 
closely related to the physical insight. Interpolation 
of current along the characteristic curves is used to 
increase accuracy. The method is relatively simple 
and can be readily programmed. The solution is 
fast and accurate and has definite advantages in 
these respects over the integro-differential equation 
approaches. 
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